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A NUMERICAL STUDY OF FORCED LITHOSPHERIC THINNING; G. Schubert, UCLA, 
C. A. Anderson, Los Alamos National Laboratory, and E. Fishbein, UCLA. 

Subsolidus lithospheric thinning by mantle plumes may be involved in the 
creation of swells, hotspots, and rifts. Among the major questions concern- 
ing this process are the timescale on which it occurs and the structure of 
the plumes. The lithosphere is known to have been substantially thinned in 
10 Ma or less' »2. 

Previous studies^ 5-6 ) 0 f the time required to thin the lithosphere by 
subsolidus convection are not in accord about the timescale of the process. 
Spohn and Schubert 5 * 6 concluded that sufficiently vigorous plumes could 
thin the lithosphere in 10 Ma, but the model only simulates the dynamics of 
the plume-lithosphere interaction in a way that does not specifically in- 
corporate the material rheology or flow pattern. Emerman and Turcotte^ 
found that the lithosphere cannot be thinned in 10 Ma, but their stagnation 
point solution cannot be matched to an external plume and boundary lqyer. 
Studies by Roberts' and Olson® indicate that within the corner region, 
at least for isoviscous vigorous convection, horizontal temperature varia- 
tions are important. For fluids with strongly temperature dependent viscos- 
ity, the overlying rigid lid should enhance these effects that are not con- 
tained in the stagnation point solutions. We have sought to clarify this 
disagreement through a numerical finite element analysis of the subsolidus 
thinning mechanism. 

Subsolidus lithospheric thinning is a process in which hot plume material 
warms the base of the lithosphere, entrains the heated lithospheric material 
in its flow, and carries the material away from the region of thinning. The 
phenomenon is basically one of forced convective heat transfer in a fluid 
with strongly temperature (and stress) dependent viscosity. Accordingly, we 
model the process as shown in Fig. 1. The parameters of the model are: T = 

temperature, u = horizontal velocity, w = vertical velocity, x = shear stress, 
a x = normal stress, t = time, Tc = surface temperature, T p = plume tempera- 
ture, w p = plume velocity, q = heat flux, and q c = conductive heat flux. At 
time t = 0 we place hot plume material with temperature T p at the base of a 
constant thickness lithospheric slab. The initial temperature in the litho- 
sphere increases linearly with depth. Plume material is forced to enter the 
base of the computational box with conditions w p , T p . The MANTLE finite 
element code^ is used to integrate the initial boundary value problem 
(heat, momentum, and continuity equations) with respect to time (and space) 
and to follow the thinning of the lithosphere (Fig. lb). Plume plus litho- 
spheric material is forced to leave the computational box along the right 
boundary as shown in Fig. lb. The viscosity of the incompressible medium is 
taken to be proportional to the exponential of its inverse absolute tempera- 
ture and to depend on the stress by a power law relation (some calculations 
were carried out with viscosity independent of stress). All other thermal, 
mechanical, and rheological properties are assumed to be constant. Litho- 
spheric thinning times will be reported as functions of rheological param- 
eters, w p , and T p . 

Morris and Canright^ 6 have argued that in steady or quasi -steady con- 
vection of a strongly temperature dependent viscosity fluid the temperature 
difference across the plume is roughly a few times the temperature difference 



134 


f ' FORCED LITHOSPHERIC THINNING 

Schubert, G. , Anderson, C., Fishbein, E. 


needed to decrease the viscosity by a factor of e (provided the rigid lid 
remains intact). This provides an important constraint on subsolidus thin- 
ning rates by convective plumes. Current studies are focused on the litho- 
spheric thinning by time-dependent plumes hypothesized by Morris to have 
large temperature differences across them. 
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Figure 1. Schematic of lithospheric thinning model. 



